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1. Introduction 

The contemporary description of the high energy evolution of QCD scattering amplitudes, 
given by the B-JIMWLK equations, has been developed over the last decade [1-20] building 
upon the pioneering ideas on gluon saturation set out by Gribov, Levin and Ryskin [21] in 
the early 80s. Notwithstanding that the B-JIMWLK set of equations has proved difficult to 
tackle and that its complete solution remains unknown, a series of numerical [22-30] and 
analytical [31-36] studies have established the asymptotic properties of the B-JIMWLK 
equations by considering their mean field limit, where the B-JIMWLK set reduces to a 
single closed equation — the Balitsky-Kovchegov (BK) equation [16-18,37,38]. Further, the 
results from this mean field equation deviate at most 10% from those obtained numerically 
from the full B-JIMWLK [30]. 

Significantly, the B-JIMWLK scheme neglects important effects [39-42] - variably 
referred to as "Pomeron loops" [42-49], "fluctuations" [41] or "wave function saturation 
effects" — due to gluon fluctuations. The insufficiencies of the B-JIMWLK framework 
become apparent once one recalls that these equations were derived under the explicit 
assumption — and, therefore, are strictly valid for such a physical situation — of large 
target gluon density and of a dilute projectile. Therefore, the B-JIMWLK formalism fails 
to properly describe the high-energy evolution for less asymmetrical systems (i.e. those 
where both target and projectile are dense). These crucial observations have led to a spurt 
of activity in this domain targeted at obtaining an evolution scheme correctly accounting 
for gluon fluctuation effects [40,42-58]. A duality transformation linking the low and high 
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density regimes [47, 48] places strict constraints on the form of any would-be complete 
evolution kernel. 

Different paths have been followed in attempting to write an evolution kernel including 
dynamics beyond B-JIMWLK. A potentially fruitful path explores the connection between 
high energy QCD evolution and a reaction-diffusion process, [41] to suggest that the full 
dynamics should be described by an equation belonging to the universality class of the 
stochastic Fisher-Kolmogorov-Petrovsky-Piscunov (sFKPP) equation [42,45,51-53,59]. A 
different strategy, the light-cone wave function approach [47-50,60], yields an evolution 
kernel not subject to the restrictions underlying the B-JIMWLK equations. This kernel, 
written in terms of the classical gluon fields of the projectile, reduces to the JIMWLK 
kernel when the fields are taken at leading order in the charge density. In this work we go 
one step further by computing the leading projectile density corrections. 

This paper is organized as follows: in Sec. 2 we briefly introduce the light-cone wave 
function approach to high energy QCD. In Sec. 3 we present our general results for the 
leading projectile density correction, and in Sec. 4 we thoroughly discuss the dipole limit 
of the results of the previous section. Our conclusions are presented in Sec. 5. 

2. Setup 

We begin by reviewing the main features of the light-cone wave function approach to the 
high energy limit of QCD as derived by Kovner, Lublinsky and Wiedemann in [50,60]. 
Consider the collision between a bunch of energetic gluons moving in the light-cone '+' 
direction — the projectile — with a dense hadronic target characterized by large gluon 
fields. Let Y be the total rapidity of the collision. The light-cone wave function of the 
incoming projectile is given by 

\yinit(Y)) = y[a?(k + ,k)}\0), (2.1) 

where a^(k + ,k) are the creation operators for gluons with color index a, longitudinal 
momentum above some hard cutoff, fc + >A + , and transverse momentum k. Assuming 
that each of the gluons in the projectile interacts independently with the target, the wave 
function after the interaction is given by 

\* f (Y)) = nS a \ Xi )a?]\0), (2.2) 

where S ab (xi) is the 5-matrix corresponding to the propagation through the target of a 
single gluon located at transverse position x%. At very high energies the interaction with 
the target eikonalizes and the S-matrix is diagonal in the transverse coordinates of the 
incoming gluons. Moreover, the S'-matrix elements depend only on the gluon fields in 
the target, A%, which allows us to take them, rather than the target fields themselves, 
as the physical degrees of freedom to describe the target. Although S ab is a quantum 
operator acting on the Hilbert space of the target fields, the commutators of 5-operators 
are suppressed by powers of the target density, which hereinafter is assumed to be very 
large. Therefore S ab can, and will, be considered as a classical c-number, in the spirit of the 
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semiclassical approach to dense gluonic systems [1-3]. With all this in mind, the scattering 
matrix of the projectile at rapidity Y can be written as 



Ey[5] = (¥(Y)\&(Y)). 



(2.3) 



High energy processes are characterized by a large separation of time scales — i.e. the time 
needed for the fast projectile to propagate through the target (the interaction time) is much 
shorter than the typical time scale under which the large, soft target fields vary. Thus, the 
fast projectile probes the target fields in a fixed, frozen configuration. Consequently, the 
physical scattering matrix is obtained by averaging over all the possible configurations of 
the gluon fields in the target: 



where the weight functional Wy[5] can be understood as the probability density for the 
target to have a certain configuration of the fields, S, at rapidity Y (see [49] for a de- 
tailed discussion of the physical meaning and properties of WyfiS]). In the wave function 
approach, the energy evolution of the system described above is achieved by boosting the 
projectile to higher rapidities, leaving the target unevolved. In this way all the information 
about the energy evolution of the system is encoded in the behavior of the projectile wave 
function as opposed to the strategy followed in the original derivation of the JIMWLK 
equation, in which the quantum fluctuations originated from boosting the target to higher 
energies were resummed in the presence of strong background fields, leading to the renor- 
malization of the weight functional, Wy . 

To first order in SY the projectile wave function at rapidity Y + 5Y is given by: 



where bi are the Weizsacker- Williams (WW) fields of the projectile, which depend uniquely 
on the projectile density operator p, defined as: 



where T a are the generators of SU(N) in the adjoint representation. For a dilute projectile 
like the one considered here, the number of gluons in its wave function is small p ~ 0(1). 
The physical meaning of Eq. (2.5) is clear: The hard, 'valence' gluons in the initial wave 
function are dressed with a cloud of soft gluons, the Weizsacker- Williams (WW) fields bi. 
These fields are determined from the classical Yang-Mills equations of motion (EOM), in 
which the hard gluons enter as an external source. The separation between soft and hard 
modes is made at an arbitrary scale A + . The boost of the projectile opens up the phase 
space for the production of new hard gluons out of the soft WW fields. This production 




(2.4) 




al a (k + ,z)}\¥(Y)), (2.5) 




(2.6) 
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process is accounted for by the last term in the rhs of Eq. (2.5). The first term corresponds 
to no gluon production, whereas the second term corresponds to virtual corrections required 
to ensure the right normalization of the wave function. 

The wave function of the projectile after the collision with the target is given by an 
analogous expression: 

\*f(Y + 6Y)) = jl- V J d 2 zb-(z,[Sp})b^z,[Sp}) (2.7) 

+ i J d 2 zb?(z,[Sp]) J* ^ ^ d ^ 1/2 S ab (z)al b (k + ,z)^ \*f(Y)), 

where now the WW fields, bi, are given by the solution of the classical EOM for rotated 
sources S ab p b . The scattering matrix of the evolved system is: 

S[5](y + SY) = (¥(Y + 6Y)\V f (Y + 5Y)). (2.8) 

From Eq. (2.3) and Eq. (2.8), it is straightforward to derive an evolution equation for 
the scattering matrix £[£] which, thanks to the Lorentz invariance of Eq. (2.4), can be 
converted into an evolution equation for the target weight functional, Wy. 

™-X*W, (2.9) 
with the kernel of the evolution given by 

X = ~hj z [pm{z ' + [Sp]mz, [Sp]) - 2b b (z, \p])b?(z, [Sp])S ba (z)]. (2.10) 

Importantly, Eq. (2.10) reduces to the JIMWLK kernel when the classical fields are 
taken at leading order in the charge density of the projectile, gp. Our goal is to derive 
higher order corrections to the JIMWLK evolution by solving the classical equations of 
motion at next order in gp. 

At this point it should be noted that the expressions in Eq. (2.5) and Eq. (2.8) for 
the evolved wave function are not complete, and they are correct only in the limit of small 
projectile density. The more general expression for the evolved wave function as given 
by Eq. (2.5) in [49] includes an extra Bogolyubov transformation with respect to the 
expression in Eqs. (2.5) and (2.8) in this paper. It is argued in [49] and assumed in this 
work that such transformation reduces to the unity operator in the limit of small projectile 
density, p — ► 0. Henceforth we restrict ourselves to the study of high density corrections 
to the kernel of the evolution arising from the expansion of the classical gluon fields, bi, in 
terms of the projectile charge density, gp. 

It is convenient to rewrite the kernel of the evolution in terms of left and right rotation 
operators, whose action on the scattering matrix of a gluonic projectile is defined as: 

J a L (x)nS] = -tr|TO(x)^y| (2.11) 
J%(x) E[S] = -tr j S(*)T a ^-y| E[5] . (2.12) 
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In terms of these operators, the kernel of the evolution reads 
3. General results 

To determine the WW fields entering the kernel of the evolution, Eq. (2.13), we need to 
solve the classical Yang-Mills equations of motion in the light-cone gauge, A + = 0, with the 
fast valence gluons playing the role of an external current. Since the hard gluons are fast 
moving in the '+' direction, the current can be written as J u (z) = g 5 u+ 5(z~ )p(z), and the 
equations read 

D lx F^(z) = r(z), (3.1) 
D^(z) = 0, (3.2) 

where Eq. (3.2) ensures the covariant conservation of the current. Since the number of 
gluons in the projectile is assumed to be small, ~ 0(1), the current that they generate is 
of order g, J ~ O(g). Thus, g will be the small parameter that controls our expansion of 
the solution. 

The general solution of these equations has been extensively discussed in [61,62]. In 
particular, it was shown in [1-3] that it is consistent to look for 'static' solutions of Eq. (3.1) 
- i.e, solutions independent of x + , with A~ = 0. Such a solution is a pure gauge with 
just transverse components. Under these assumptions Eq. (3.2) is trivially satisfied and, by 
writing the transverse components of the field as Ai{z~ , z) = 0(z~)bi(z), Eq. (3.1) becomes 

dib?(z) +gf ahc b b l {z)b<i{z) = gp a (z) , (3.3) 
diVfrz) - dtfW+gf^bfcWjiz) = 0, (3.4) 

Expanding the solutions in powers of g: 

b?(z)=g(b? 1 (z)+g 2 b° 2 (z) + 0(g i )), (3.5) 

we get 

b * i{z) = hi d2x ^w pa{z) = hl Sx & x )p a ^ < 3 - 6 ) 

= -J^f^j d 2 xd 2 y(Xd t Y)p\x)p c (y), (3.7) 

where we have made use of the following shorthand notation for the coordinate dependence 
of the solution: 

diX = dfln(\z-x\X) , 
XdiY = X(dfY)-(dfX)Y, (3.8) 

where df denotes the partial derivative with respect to the transverse components of z, 
i = l,2. 
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The expansion of the solution in Eq. (3.5) immediately translates into an expansion 
for the kernel of the evolution in powers of a s = g 2 /Air: 

X = a s ( X 1) + u sX 2) + 0{a 2 s )), (3.9) 

whose leading term is the JIMWLK kernel: 

asX l) = JIMWLK = JjdaMY) [J a L (z)Jl(z) 

+ r R (z)J a R {z) - 2Jt(z)J b R (z)S ab (z)] . (3.10) 

The first correction to JIMWLK, \ 2 \ is obtained by keeping terms up to order a 2 s in the 
product WW fields, yielding 

X 2) = "7^2 / f abc (0^)(Xd t W) \jt(x)jt(y)Ji{w) + JUx)J R (y)J R (w) 

^r^XXYdtW^iz) \j b L (x)J R (y)J d R (w) + J b L (x)J c L (y)J d R W] ■ (3-11) 



Note that from the second order solution of the EOM we could immediately derive 
part of the 0(a^ ) corrections to the kernel. However, a complete derivation of the 0{o? s ) 
term in Eq. (3.9) would require the 0(g 5 ) solution to the EOM, which is beyond the scope 
of this paper. 

4. Dipole model limit 

The color structure of a generic projectile composed of gluons can be rather complicated, 
consisting of different color multipoles mix with each other through the evolution in a 
highly non-trivial way. However, in the large- N limit this complicated structure is greatly 
simplified and the high energy evolution can be recast in terms of dipole degrees of freedom. 
More precisely, the JIMWLK equation is equivalent to an infinite hierarchy of coupled 
differential equations for the correlators of the gluon fields. In the large- N limit the whole 
hierarchy decouples and one is left with a single, closed, non-linear evolution equation for 
the dipole scattering amplitude — the BK equation [16-18,37,38]. 

In this section we explore the color structure of the correction to JIMWLK evolution 
derived in the previous section. In order to do so, we consider an initial projectile entirely 
describable by dipole degrees of freedom: 

E = £[*], (4.1) 

where 

s(x,y) = -tv{sUy)S F (x)} (4.2) 

is the scattering matrix for a q-q dipole, with x and y the transverse coordinates of, re- 
spectively, the quark and the antiquark. The subscript F in the rhs of Eq. (4.2) indicates 
that the scattering matrix is to be taken for particles in the fundamental representation of 
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SU(N), and will be usually omitted in the following. Our purpose is to study the resulting 
color structure of the projectile under evolution, i.e. we want to calculate x ^[ s ]- The 
action of the left and right rotation operators on a dipole-like projectile scattering matrix 
is given by: 



J a L (x)s[u,v] = -I t r{«(x)^^ 
J%(x)s[u,v] = -l tT { S ( x ) ta js^) 



S^(x)t a 
t a S\x) 



SS*{x) 
5 

5S*(x) 



s[u,v] , 
s[u,v] , 



(4.3) 
(4.4) 



where t a are now the SU(N) generators in the fundamental representation. From Eqs. 
(4.3) and (4.4) it can be proved that both Jl and Jr are hermitian operators, a property 
which immediately ensures the hermiticity of the kernel of the evolution. 

4.1 General expressions 

Using the following notation for the functional derivatives of the scattering matrix: 

6 5 



Juv;...;qt 



s = 



Ss(u,v) 5s(q,t) 



S[ S ], 



(4.5) 



and defining S z = Sf(z), we can write the action of the different pieces of the kernel in 
Eq. (3.11) on E[a] as 

r bc J a R (x)J R (y)J R (w)E[s] = -lf^[S(v-d)-S^-w)]{ (4.6) 
[S(v - y)+5{u-y)] [S(v -x)- 5{u - x)]tr {Sj S v t a t b t c }5 uv Z 
+ L[5(v-y)+5(u-y)][5(r-x)-5(p-x)]tr{SiSjt c }tr{S^ 

+^[S(r-y)-S^-y)][S(v-x)+8(u-x)]ti{SlS v tH c }^ 

+^[<5(r-y)-<5fc-y)][<5(r-x)+<5(^ 

+^[S(r-y)-S^-y)}[S(t-x)-S(q-x)}ti{SiS v t c }tT{S^ 



f abc J a L (x)J b (y)J c L (w)X[s] = -^f abc [5{v-w)-5{u-w)]{ (4.7) 
[S(v - y)+S(u - y)} [S(v -x)- 5{u - x)]tr {Stt c t b t a S V }5 UV Z 
+ ^[5(«-y)+5(n-y)][<5(r-x)-^-x)]tr{5t t ^4tr{S'tr5JV;™S 

+ ^[5(^-y)-^-2/)]^(^-^)+^(^-^)]tr{^ a 54tr{4t 6 5 r }<5 pr ;«^ 

+ ^[5(r-y)-5fr-y)][5(r-x)+5(p-x)]tr{SifS v }tr{S^ 

+ ^[S(r-y)-5^-y)][5(t-x)-5(q-x)]tr{Slt c S v }tr{Sy 



-7- 



r d J b L {x)J c R {y)J d R {w)ns] = -±r d [5(v-w)-5(u-w)}{ (4.8) 
[5(v-y)+5(u-y)][5(v-x)-5(u-x)]tv{Slt b S v t c t d }5 uv J: 

+^[*h!/)-^-!/)][i(^-x)4( U -.)]tr{5tt t y>{5t5 r f}V;»S 
+^l5(r-y)-5fa-y)}[5(r-x)-5(p-x)}ti{SiS v t d }tr{Slt b ^ 

+^(r-y)-*fr-y)M*-s)-<^ 

and 

f hcd Ji(w)J d L (y)J a R {x)ns] = -±f bdc [S(v-x)-5(u-x)]{ (4.9) 
[5(v-w)-5(u-w)}[5(v-y)+5(u-y)]tr{Sit d t c S v t a }5 uv Z 
+^[5(v-w)-5(u-w)][5(r-y)-5^-y)]tT{Slt d S v t a }tT{S^ 



+l[$(r-t U )- ( Sfr- w )P(t,-^-<^ 

+4[^-^)-^(p-^][5(r-y)+<5(p-y)]tr{5t5X}tr{4^ c 5J5 pr ;™S 



+^(r-™)-<5(p-™)P(t-y)-<5( (7 ^ 



4.2 Diagrammatic interpretation 

Despite the fact that the expressions derived in the previous section appear complicated, 
they allow for a very clear physical interpretation in terms of diagrams according to the 
following rules: the action of the left (right) rotation operator, J^my 011 the projectile 
scattering matrix brings in a new dipole, along with the the corresponding 1/N suppression 
factor, which emits a new gluon of color a before (or after) the interaction with the target. 
Such emission may happen either from the quark line or from the antiquark. In the former 
case a, relative minus sign is picked up. Subsequent actions of the rotation operators may 
act either on bringing new dipoles into the diagram, or on the preexisting dipole, 

which emits a new gluon. These rules are sketched in Fig. 1. 

In order to fully describe the diagrams in terms of fundamental constituents, quark 
and antiquark lines, we make use of the Fierz identity: 

= \(Su5jk ~ jfSijhi), (4.10) 

which can be translated into diagrammatic language substituting the gluon lines by quark- 
antiquark lines as indicated in Fig. 2. 

The propagation of each quark (antiquark) line through the target at transverse co- 
ordinate x is accounted for the scattering matrix, Sf(x) (S f (x)). Finally, the trace over 
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Figure 1: Diagrammatic representation of J£E (top) and J£J£,£ (bottom). The target is repre- 
sented by the vertical thick line. 
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Figure 2: Pictorial representation of the Fierz identity Eq. (4.10). 



closed fermion lines has to be taken. Following these rules, it is straightforward to rederive 
the results in Eqs. (4.6)-(4.9). As an example, we show in Fig. 3 the diagrams correspond- 
ing to the action of JlJlJl^[s], Eq. (4.7), and JlJlJr^[s], Eq. (4.9), in which, for the 
sake of simplicity, we have kept the gluon lines in the adjoint representation without using 
the Fierz identity. 

4.3 Symmetry considerations 

Before proceeding further it is convenient to divide our final result into three pieces ac- 
cording to their leading power of 1/N as given by the number of dipoles appearing in the 
relevant diagram, see the previous subsection: 

x 2) s[ s ] = (xf /N + xf /N 2 + xf /N3 ) £[4 (4.11) 

To complete the calculation we still have to plug in the results from Eqs. (4.6)-(4.9). 
into Eq. (3.11), perform the contraction of the color indices and integrate over x,y,w. 
Further, the cross terms LRR, Eq. (4.8), and LLR, Eq. (4.9), must be contracted with the 
scattering matrix for a single gluon, S^, which can be rewritten in terms of matrices in 
the fundamental representation by means of the Fierz identity: 

Sf(z) = 2tr{S^ F (z)t b S F (z)t a }. (4.12) 

The SU (N) algebra can be worked out using the relations listed in the Appendix. 

After all this is done, we note that the contribution to the kernel coming from the 
pieces of order 1/N in Eqs. (4.6) and (4.7) cancel each other out. An analogous cancellation 
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Figure 3: Diagrams corresponding to the LLL (plots on the left) and LLR (plots on the right), see 
the accompanying text. The diagrams include one (ex 1/N, top), two (oc 1/N 2 , middle) and three 
(tx 1/iV 3 , bottom) active dipoles. The thick vertical line represents the target. 



happens between the leading 1/N contributions from Eqs. (4.8) and (4.9). Therefore 

X? /JV S[ S ] = 0. (4.13) 

It is clear from the diagrammatic rules derived in 4.2 that this cancellation corresponds 
to the diagrams in which all the emissions and absorptions of gluons happen in a single 
dipole, while the other dipoles in the wave function remain as spectators. More precisely, 
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the top diagrams in Fig. 3 cancel against their respective complex conjugate diagrams (not 
shown in the figure). 

An analogous cancellation occurs for the 1/N 2 terms in Eqs. (4.6) and (4.7). Dia- 
grammatically, the middle-left diagram in Fig. 3 cancels against its complex conjugate. 
Therefore, the leading 1/N 2 contribution is given by the 1/N 2 terms in Eqs. (4.8) and 
(4.9) (middle-right diagram in Fig. 3 plus its complex conjugate). We get 



xf /N2 m- (^)2^2 j/ 



d 2 z 



- N[d t (R-P)][U S V][tr{StS v StS r SlS z } - -tr{SlS v }tT{Sj,S r }] 



N 

- X - [d t (V-U)][(R-P) S (F-^lttrjSt^^M.St^l-trl^t^^t^M^t^}] 
+ l - [di{R-P)\[(R-P) di (V-U)][ti{SiS r SlS v }tr{SlS z } - tr{SlS z SlS v }tr{SiS r }] 

- \ [di(V-U)][(R-P) S (F-^lttrlS^t^^M^^J-trl^^Sl^MSt^}] 
+ \ mV-U)][{V-U) di (R-P)][ti{S z StS r Sl}ti{SiS v } - tv{sls v sis r }tv{sls z }\ 

+ N[di{V-U)][P di R][t?{SiS v S\S r SlS z } - ^{SlS v }tr{SlS r }] ^ 



6s(p, r)5s(u, v) 

(4.14) 



This result can be further simplified by noting that any wave function or weight func- 
tional of a gluonic/dipole configuration has to be completely symmetric under the exchange 
of any number of gluons/dipoles, such that the exchange (u,v) <-> (p,r) leaves the action 
of the functional derivative S uv;pr T, unchanged. Under such an exchange, many terms in 
Eq. (4.14) cancel each other, yielding: 



[diiR-P^UdiV) - d(V-U)(PdiR)][Ntr{StS v SiS r SlS z }] 
+d l (V-U)(R-P) di (V-U) \tr{StS v SiS r }tr{SlS z }-tr{SiS z SlS r }tr{SiS^ 



+tv{SlS z StS r }ti{SiS v } - tv{SlS v SiS r }tv{SlS z }] } — ^ 



Ss(p, r)5s(u, v) 



(4.15) 



Calculating the 1/iV 3 terms in an analogous way we get: 
i 1 1 



{d i (R-P)}{(T-Q)d l (V-U)} 



2 (2vr) 2 iV 3 

^{S^S v SqStSpS r } + ^{S^SfS^SrSpSy} + ti{Sj,S r SpS z Sl l S v SqSt} 

<5 3 £[s] 



+tv{sls r sts v sis t sls z } 



5s(q, t)5s(p, r)8s(u, v) 



(4.16) 
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Eqs. (4.13), (4.15) and (4.16) are the central result of this paper. One of their main 
features is the fact that the action of the kernel on a dipole-like projectile cannot be entirely 
recast in terms of dipole degrees of freedom. On the contrary, the evolution generates a 
complicated color structure consisting in the mixing of dipoles, quadrupoles, sextupoles 
and octupoles as given by these equations, which can only be expressed partially as dipole 
degrees of freedom by increasing the power in \/N . Consequently, these corrections to the 
leading JIMWLK kernel bring no corrections to the mean field BK equation. 



5. Conclusions 

In this work we have calculated the 0(at 2 s ) corrections to the JIMWLK kernel that arise from 
the 0(g 3 ) solutions to the classical Yang-Mills equations. These are the first corrections in 
the charge density of the projectile to JIMWLK evolution, and therefore partially account 
for the coherence effects in the projectile gluon emission which drives small- x evolution. In 
the context of the present discussions on the quest for evolution equations for scattering of 
a dense projectile which would contain the so-called pomeron loops, our result accounts for 
the part of these corrections not arising from the non-commutativity of the target fields. 
Thus, they are restricted to the dilute-dense scattering situation. The same systematic 
technique could be used to improve the results, although admittedly the computation of 
even higher orders looks technically challenging. 

Our main result is the cancellation of the leading 1/N contributions, together with 
Eqs. (4.15)-(4.16). We provide a diagrammatic interpretation of these results. These 
corrections are sub leading in 1/N and exhibit a complicated color structure. They do not, 
therefore, provide any correction to the BK equation but, rather, add new terms to the 
Balitsky hierarchy. 
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A. Some color algebra 

In this appendix we list some relations needed to derive the results in section 4. Using the 
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basic SU(N) relations 

jade jbef jcfd _ jabc ^ 
d ade d bef f cfd = ^_ _ ^ jabc ^ 

d adejbeffcfd = __ jabc ^ ( A 3 ) 

d ade d bef d cfd = ^ jabc ^ ^ ^ 

tv{t a t b t c t d } = -i- <5 ab <f d + ^ , (A.5) 

47V 8 

where we have introduced the notation: 

g abc = 4 tr { t ^fe t c| = ^abc + if abc^ . -abc _ ^abc _ -ya&c) _ ^ g) 

Since the generators of SU (TV), t a , along with the unit matrix form a basis of the matrices in 
the fundamental representation, we expand any arbitrary matrices, X,Y . . . in the following 

way 

X = x + x a 1 t a , Y = y + y1t\... (A.7) 
With the help of the relations listed above, we get 

f abc tr{XYt a t b t c } = j(N 2 - l)tr{XY}, (A.8) 

f ahc tr{Xt a Yt b t c } = ^ (tr{X}tr{Y} - ^tr{XY}^j , (A.9) 

f abc tv { Xt a Yt b }tr{Zt c } = l - [tr{XZ}tr{Y} - tr{YZ}tr{X}] , (A.10) 



f aoc tr{XYt a f}tr{Zt c } = — [tr{XY Z] - —tr{XY}tr{Z} j , (A.ll) 

f abc tr{Xt a }tr{Yt b }tr{Zt c } = % - [tr{X ZY} - tr{XYZ}] . (A.12) 

These are the relations required to perform the contraction of the color indices in Eqs. 
(4.6)-(4.9). 
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